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ABSTRACT 

Using analytical methods and Monte Carlo simulations, we analyze a new 
statistic designed to detect isolated step-like discontinuities which are coherent 
over large areas of Cosmic Microwave Background (CMB) pixel maps. Such 
coherent temperature discontinuities are predicted by the Kaiser-Stebbins 
effect to form due to long cosmic strings present in our present horizon. The 
background of the coherent step-like seed is assumed to be a scale invariant 
Gaussian random field which could have been produced by a superposition 
of seeds on smaller scales and/or by inflationary quantum fluctuations. The 
effects of uncorrelated Gaussian random noise are also considered. The 
statistical variable considered is the Sample Mean Difference (SMD) between 
large neighbouring sectors of CMB maps, separated by a straight line in two 
dimensional maps and a point in one dimensional maps. We find that including 
noise, the SMD statistics can detect at the la to 2a level the presense of a 
long string with GiJ,{vsjs) — ^(y'^^'^^ — ^-^ ^ while more conventional 
statistics like the skewness or the kurtosis require a value of Gn almost an order 
of magnitude larger for detectability at a comparable level. 

Subject headings: cosmic microwave background - cosmic strings 
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1. Introduction 

The major progress achieved during the past 15 years in both theory and cosmological 
observations has turned the search for the origin of cosmic structure into one of the 
most exciting fields of scientific research (for a good review see Efstathiou 1989). Despite 
the severe constraints imposed by detailed observational data on theories for structure 
formation the central question remains open: What is the origin of primordial fluctuations 
that gave rise to structure in the universe? Two classes of theories attempting to answer 
this question have emerged during the past ten years and have managed to survive through 
the observational constraints with only minor adjustments. 

According to the first class, primordial fiuctuations are produced by quantum 
fluctuations of a hnearly coupled scalar fleld during a period of inflation (Hawking 1982; 
Starobinsky 1982; Guth & Pi 1982; Bardeen, Steinhardt & Turner 1983). These fluctuations 
are subsequently expected to become classical and provide the progenitors of structure in 
the universe. Because of the extremely small linear couphng of the scalar fleld, needed to 
preserve the observed large scale homogeneity, the inflationary perturbations are expected 
by the central limit theorem, to obey Gaussian statistics. This is not the case for the second 
class of theories. 

According to the second class of theories (Kibble 1976; Vilenkin 1981; Vilenkin 1985; 
Turok 1989; Brandenberger 1992; Perivolaropoulos 1994), primordial perturbations are 
provided by seeds of trapped energy density produced during symmetry breaking phase 
transitions in the early universe. Such symmetry breaking is predicted by Grand Unifled 
Theories (GUT's) to occur at early times as the universe cools and expands. The geometry 
of the produced seeds, known as topological defects is determined by the topology of the 
vaccuum manifold of the physically realized GUT. Thus the defects may be pointhke 
(monopoles), linelike (cosmic strings), planar (domain walls) or collapsing pointlike 
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(textures) . 

The cosmic string theory (Vilenkin 1981) for structure formation is the oldest and 
(together with textures (Turok 1989)) best studied theory of the topological defect class. By 
fixing its single free parameter Gji (// is the effective mass per unit length of the wiggly string 
and G is Newtons constant) to a value consistent with microphysical requirements coming 
from GUT's, the theory may automatically account for large scale filaments and sheets 
(Vachaspati 1986; Stebbins et. al. 1987; Perivolaropoulos, Brandenberger & Stebbins 1990; 
Vachaspati & Vilenkin 1991; Vollick 1992; Hara & Miyoshi 1993), galaxy formation at epochs 
2; ~ 2 — 3 (Brandenberger et. al. 1987) and galactic magnetic fields (Vachaspati 1992b). 
It can also provide large scale pecuhar velocities (Vachaspati 1992a; Perivolaropoulos & 
Vachaspati 1994) and is consistent with the amplitude, spectral index (Bouchet, Bennett 
& Stebbins 1988; Bennett, Stebbins & Bouchet 1992; Perivolaropoulos 1993a) and the 
statistics (Gott et. al. 1990; Perivolaropoulos 1993b; Moessner, Perivolaropoulos & 
Brandenberger 1994; Coulson et. al. 1994; Luo 1994; Magueijo 1995b) of the cosmic 
microwave background (CMB) anisotropics measured by the COBE collaboration (Smoot 
et. al. 1992; Wright et. al. 1992) on large angular scales {9 - 10°). Other planned CMB 
experiments (see e.g. MAP 1997, COBRAS/SAMBA 1997, see also the review by Scott et. 
al. 1995) of equally high quality but on smaller angular scales are expected to provide a 
wealth of information within the next few years. 

The CMB observations provide a valuable direct probe for identifying signatures of 
cosmic strings. The main mechanism by which strings can produce CMB fluctuations on 
angular scales larger than 1-2 degrees has been well studied both analytically (Brandenberger 
& Turok 1986; Stebbins 1988; Veeraraghavan & Stebbins 1990; Perivolaropoulos 1993a; 
Perivolaropoulos 1993b; Moessner et. al. 1994) and using numerical simulations (Bouchet 
et.al. 1988; Bennett et. al. 1992) and is known as the Kaiser- Stebbins effect (Kaiser & 
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Stebbins 1984; Gott 1985). According to this effect, moving long strings present between 
tfie time of recombination tree and the present time to, produce step-hke temperature 
discontinuities between photons that reach the observer through opposite sides of the string. 
These discontinuities are due to the pecuhar nature of the spacetime around a long string 
which even though is locally fiat, globally has the geometry of a cone with deficit angle 
SttG/x. The magnitude of the discontinuity is proportional to the deficit angle, to the string 
velocity Vg and depends on the relative orientation between the unit vector along the string 
s and the unit photon wave-vector k. It is given by (Stebbins 1988) 

ST 

— = ±4:nGfxVs-fsk ■ {vs X s) (1) 

where 7^ is the relativistic Lorentz factor and the sign changes when the string is crossed. 
The angular scale over which this discontinuity persists is given by the radius of curvature 
of the string which according to simulations (Bennett & Bouchet 1988; Allen & Shellard 
1990; Albrecht & Stebbins 1993) is approximately equal to the horizon scale. The growth 
of the horizon from tree to to results in a superposition of a large number of step-like 
temperature seeds of all sizes starting from about 2° (the angular size of the horizon at tree) 
to about 180° (the present horizon scale). By the central limit theorem this large number 
of superposed seeds results in a pattern of fluctuations that obeys Gaussian statistics. 
Thus the probabihty distribution for the temperature of each pixel of a CMB map with 
resolution larger than about l" — 2" is a Gaussian (Allen et. al. 1996, Coulson et. al. 
1994, Perivolaropoulos 1993b; Perivolaropoulos 1993c). It has therefore been considered 
to be impossible to distinguish structure formation models based on cosmic strings from 
corresponding models based on inflation, using CMB maps with resolution angle larger 
than 1° — 2'^ (Ferreira & Magueijo 1997) . Theoretical studies have therefore focused on 
identifying the statistical signatures of cosmic strings on angular scales less than 1° (Turok 
1996, Coulson et. al. 1994) where the number of superposed seeds is smaller and therefore 
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the non-Gaussian character of fluctuations is expected to be stronger [] 

These efforts however have been faced with the comphcated and model dependent 
physical processes occuring on small angular scales. Such effects include isolated foreground 
point sources, recombination physics, string properties on small scales (kinks, loops etc) 
which require detailed simulations of both the string network and the cosmic background, in 
order to be properly taken into account. Even though preliminary efforts for such detailed 
simulations are in progress (Allen et. al. 1996), it has become clear that it will take some 
time before theory and experiments on angular scales less than a few arcmin reach accuracy 
levels leading to detectable non-Gaussian string signatures. 

An alternative approach to the problem is instead of focusing on small scales where 
the number of superposed seeds is small, to focus on larger angular scales where despite 
the large number of superposed seeds there is also coherence of induced fluctuations on 
large angular scales. Fluctuations on these scales may be viewed as a superposition of a 
Gaussian scale invariant background coming mainly from small scale seeds plus a small 
number of step-like discontinuities which are coherent and persist on angular scales larger 
than 100*^. These are produced by long strings present in our present horizon. Our goal 
is to find a statistic optimized to detect this large scale coherence and use it to find the 
minimum amplitude of a step function that can be detected at the la level relatively to 
a given scale invariant Gaussian background. Such a statistic is equally effective on any 
angular resolution scale and its effectiveness is only diminished as the number of pixels of 
the CMB map is reduced or the noise is increased. The statistical variable we focus on, in 
what follows is the Sample Mean Difference (SMD) between large neighbouring sectors of a 

^The non-Gaussian features for texture maps are stronger than thosed of cosmic strings 
mainly because of the generically smaller number of textures per horizon volume (Gangui 
1996; Magueijo 1995a). 
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CMB map. These sectors are separated by a random straight hne in two dimensional maps 
or by a random point in one dimensional maps. The union of the two sectors gives back the 
complete map. We show that the statistics of the SMD variable are much more sensitive 
in detecting the presence of a coherent step-like seed than conventional statistics like the 
skewness or the kurtosis. 

The structure of this paper is the following: In the next section titled 'Sample Mean 
Difference' we study analytically the statistics of the SMD variable and show that its 
average value is a sensitive quantity in detecting the presence of a randomly positioned 
step-function on top of a gaussian map. We then compare with the sensitivity of the 
statistics skewness and kurtosis. We find that the sensitivity of the SMD statistics is 
significantly superior to that of skewness and kurtosis in detecting the step function. These 
analytical results are shown for the case of one-dimensional maps but the extension to the 
case of two dimensional maps is straightforward. 

In the third section titled 'Monte Carlo Simulations' we perform Monte Carlo 
simulations of Gaussian maps with flat and scale invariant power spectra, with and without 
step-like discontinuities in one and two dimensions. We also include noise with signal 
to noise ration ^ = 1.0. Applying the statistics skewness, kurtosis and average of SMD 
on these maps we verify the results of section 2 and find the minimum step-function 
amphtude that is detectable by the average SMD statistic. Finally is section 4 we conclude, 
summarise and discuss the prospect of applying the mean of SMD statistic to presently 
available CMB maps including the COBE results. That analysis is currently in progress 
(Athanasiou,Perivolaropoulos & Simatos 1997). 
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2. Sample Mean Difference 

Consider an one dimensional array of n pixel variables Xn. Let these variables be 
initially distributed according to a standardised Gaussian probability distribution. Consider 
now a step-function of amplitude 2a superposed so that the discontinuity is between pixels 
io and io + 1 (Fig. 1). The new probability distribution for a random pixel variable x is 

P(x) = ^=e — + —r=^e — (2 



where f = ^- We are looking for a statistic that will optimally distinguish between a 
Gaussian array with a superposed step-function and a Gaussian array without one. The 
obvious statistics to try first are the moments of the distribution (2) with a — and a ^ 0. 

The moment generating function corresponding to (2) is: 

M{t) = /e"*+^ + (1 - /)e-"*+^ (3) 

The mean /^(q;,/), variance a^{a,f), skewness s{a, f) and kurtosis k{a, f) can be obtained 
in a straightforward way by proper differentiation of M{t) as follows: 

//(a,/) = <X>^af-a{l-f) (4) 
a'iaj) = <{x-pif>^l + Wf{l-f) (5) 

a 

3 + 8aV(3 + 2^2 - 3/2 - 8«V + 12^2/^ - ^c^^f) 



(l + 4aV(l-/)r 

For q; = we obtain the Gaussian values for the skewness and the kurtosis s(0, /) = 
, A;(0, /) = 3 as expected. For a ^ the moments deviate from the Gaussian values. In 
order to find the minimum value of a for which the moments can distinguish between a 
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Gaussian pattern and a Gaussian+Step pattern we must compare the deviation of moments 
from their Gaussian values with the standard deviation of the sample moments. The mean 
values of the skewness and the kurtosis are easily obtained by integrating with respect to / 
from to 1 i.e. assuming that it is equally probable for the step-function to be superposed 
at any point of the lattice. 

s{a) = < s{a, /) >= / dfs{a, /) = (7) 

k{a) = <k{a,f)>= f'dfkiaj) (8) 

Jo 

These values are to be compared with the standard deviations of the moments, obtained as 
follows: The variance of the skewness over several n-pixel array realizations with fixed / 
and a is 

As'{aJ)^<{s-sf> (9) 

where s = is the sample skewness from a given pixel array realization, s is the 

actual skewness and Si = i^^i^^. Now 

n< si> 

< s >= si >= s (10) 

n 



Also 



< >= - < s2 > +(1 - -) < Sj >2 (11) 



where j any pixel number (j e [l,?^])- Thus 



As\a, /) = 1(< ^2 > _ < >2) - -4 < i^J - f^T > (12) 



Similarly for the variance of the sample kurtosis we have 

Ae{a, /) = 1(< fcj > - < kj >') (13) 

with kj — -^{xj — < ^"j >— ^ < {^3 ~ A*)^ >■ It is straightforward to obtain all the 

above moments by differentiating the generating functional and using 
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Now the minimum value amin of a detectable at la level is obtained from the equations 

/ df[s{amin, f) - ^S{amvn, /)] = (15) 

Jo 

f'df[{k{armnJ)-^)- thiamin J)] = (16) 

Jo 

Since (from eq. (7)) s{a) — which is equal to the Gaussian value, the skewness can 
only be used to detect a step function by comparing the standard deviation As for a = 
and a ^ 0. By demanding As{amin) < 2As{a = 0) we obtain amin < 2.5. This result 
is independent of the number of pixels n. For the kurtosis we obtain from eqs. (13, 16) 
amin — 4 for n = 10^ while for amin = 0.5, n ~ 10^ is required. 

Using the alternative test i.e. demanding Ak{ajnin) > 2Ak{a — 0) we obtain a > 2 
and this result is independent of the number of pixels n as in the case of skewness. Thus 
for the usual pixel maps where n is up to 0(1000) the kurtosis is not able to detect a step 
function with a < 2 at the la level. As in all cases discussed in this paper a is measured 
in units of standard deviation (rms) of the underlying Gaussian map. This result remains 
unchanged for other statistical variables defined by local linear combinations of pixels (e.g. 
differences of neighbouring pixel variables (Moessner et. al. 1994, Coulson et. al. 1994)) 
since the effect of a single discontinuity remains negligible if the long range coherence is not 
taken into account. 

For CMB temperature maps with (^)rms — 2 x 10~^ the detectable value of Gfi is 

a = 4:TTGii{v,-fs) cos ^ > 4 X 10"^ =^ Gii{vs-fs) cos ^ > 4 x 10"^ (17) 

where 6 is an angle obtained from the relative orientation of the string with respect to the 
observer. According to simulations < Vsjs >rms— 0-2 and for < 2 x 10~^ the detection 
of the Kaiser-Stebbins effect using statistics based on skewness and kurtosis is not possible. 
This excluded range however includes all the cosmologically interesting values of G^. 

It is therefore important to look for alternative statistical variables that are more 
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sensitive in detecting the presence of coherent discontinuities superposed on Gaussian maps. 
As we will show, the Sample Mean Difference (SMD) is such a statistical variable. 

Consider a pixel array (Fig. 1) of n pixel Gaussian random variables Xj with a step 
function covering the whole array, superposed such that the discontinuity is located just 
after pixel io. To every pixel k of the array we may associate the random variable defined 
as the difference between the mean value of the pixels 1 through k minus the mean value of 
the pixels k + 1 through n. It is straightforward to show that 

Yk = AXk + 2a'^^—^ ke[l,io] (18) 
n — K 

Yk = AXk + ke [to, n-1] (19) 

where AXk = jY,']=iXj — :^jyj=k_^.iXj. Thus we have constructed a new array Y^, 
{k = 1, ...,n — 1) from the sample mean differences (SMD) of the original array. We will 
focus on the average value Z of the SMD defined as: 

n— 1 



-rE'^k (20) 

k=l 

Using eqs. (19,20) we obtain 

With the definitions / = io/n and ^ = k/n and the assumption n >> 1 we obtain: 

Z = J2 - 2a[(l - /) ln(l - /) + / In /] (22) 

k=i 

Thus the mean of Z over many realizations of the array is 

< Z >= J2 < > -H t dfflnf]^a (23) 

The variance of Z is due both to the underlying Gaussian map and to the variation of 
f — io/n (assuming a fixed). The variance due to the gaussian background is 
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where e — 0{^), ^ — k/n, n » 1 and we have used the fact that the variance of the sample 
mean of a standardized Gaussian population with size J is j- Now from eq. (24) we obtain 



9 , 9 21nn 

(7i%~-elne2~— — (25) 

The variance of the /-dependent part of Z is 

(T2,Z =< ^2 > - < Z2 >^ (26) 

where Z2 = -2a[(l - /) ln(l - /) + / In /]. From eq. (23) we have < Z2 >= a and < Z| > 
is easily obtained as 

<Zi>^ £dfZl(f)^^a (27) 

Thus 

4 = ^Iz + ^Iz = ^ + (28) 
In order to be able to distinguish between a Gaussian+Step map and a purely Gaussian 
one, at the Icr level we demand that 

< Z >„^o - <Z >a=o> (yz (29) 
This implies that the minimum value of a, otmin that can be detected using this test is 

oim^ = (^)V2 (30) 



n 

and for n = 0(10"^) we obtain amin — 0.2 which is about an order of magnitude improvement 
over the corresponding sensitivity of tests based on the moments skewness and kurtosis. 
The reason for this significant improvement is the fact that the SMD statistical variable 
picks up the coherence properties introduced by the step function on the Gaussian map. 
The moments on the other hand pick up only local properties of the pixels and do not 
amplify the long range coherence of the step-like discontinuity. 

Our analysis so far has assumed that the Gaussian variables Xj are independent and 
that the only correlation is introduced by the superposed step-function. In a realistic 
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setup however the underlying Gaussian map will be scale invariant and thus there will be 
correlations among the pixels. These correlations will also be affected by the instrument 
noise. In addition, our analysis has been limited so far to one dimensional maps while 
most CMB experiments are now obtaining two-dimensional maps. In order to take all 
these effects into account we need to apply the statistics of the SMD variable onto maps 
constructed by Monte Carlo simulations. This is the focus of the following section. 



3. Monte-Carlo Simulations 

We start by constructing an array of n Gaussian random variables Xj, j = 1, ...,n with 
a power spectrum P{k) — k~'^. Thus the values Xj associated with the pixel j is obtained 
as the Fourier transform of a function g{k) {k — 1, ...,n) with the following properties: 

• For each fc, the amplitude \g{k)\ is an independent random variable with mean and 
variance P{k) = 1/k"^. 

• The phase 6k of each Fourier component g{k) is an independednt random variable in 
the range [0, 27r] with uniform probability distribution P{9k) — 

• The Fourier components are related by complex conjugation relations neeeded to give 
a rea/ variable Xj. 

The discrete Fourier transform definition used is 

k=l 

and the numerical programming was implemented using Mathematica (Wolfram 1991). In 
order to have real Xj, the conditions 

Img(O) = (32) 
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g{k + 2) = g*{n-k) k^0,...,n-2 (33) 
Img(^ + 1) = (34) 

must be satisfied. The array Xj obtained in the way described above is then standardized 
to a new array X^, with 

X; . (36) 

where /j, and are the sample mean and sample variance for the realization of the array 
Xj. A new array X^ is then constructed by superposing to the array Xj a step function of 
amplitude 2a with discontinuity at a random point iq. The array X^ is thus obtained as 

X;=X; + alr^, j = l,...,n (36) 

^ \J-^o\ 

Next we apply the statistics discussed in the previous section to several realizations of the 
arrays Xj and X' in an cfTort to find the most sensitive statistic that can distinguish among 
them. Our goal is to also find the minimum value of a that can be distinguished by that 
statistic at the la level, thus testing the analytical results of the previous section. 

We have used a lattice with 2000 pixels and a scale invariant power spectrum which 
for one-dimesional data is P{k) — k~^. In Table 1 we show the results for the skewness, 
the kurtosis and the average SMD for the Xj arrays, with a =0, 0.25, 0.50 and 1.0. The 
SMD average was obtained as in section 2 by first constructing the array of sample mean 
differences and then obtaining its average value, predicted to be equal to a by the analytical 
study of section 2. 

These statistics were applied to 50 random realizations of the array X?. The mean 
values of the statistics considered with their la standard deviations obtained over these 50 
realizations are shown in the following Table 1. 
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Table 1: A comparison of the effectiveness of the statistics considered, in detecting 
the presence of a coherent step discontinuity with amphtude 2q; relative to the standard 
deviation of the underlying Gaussian map. No noise was included in these simulations and 
the full map was used in obtaining the SMD average. 



a 


Skewness 


Kurtosis 


SMD Average 


0.00 


0.01 ±0.11 


2.97 ±0.19 


0.02 ±0.31 


0.25 


0.01 ±0.11 


2.95 ±0.20 


0.25 ±0.33 


0.50 


0.02 ±0.11 


2.88 ±0.21 


0.48 ±0.38 


1.00 


0.03 ±0.20 


2.82 ± 0.32 


0.98 ± 0.48 



The analytical prediction of section 2 for the SMD average value a is in good agreement 
with the results of the Monte Carlo simulations. The standard deviation of this result is 
not in such a good agreemnent with the analytical prediction because the assumption of 
complete independence among pixels made by the analytical treatment is not realized in the 
Monte Carlo simulations where a scale invariant spectrum was considered and thus there 
was a non-trivial correlation among the pixels of the arrays. 

The effects of adding uncorrelated Gaussian noise with signal to noise ratio ^ = 1 are 
shown in Table 2. This table was constructed by adding an uncorrelated Gaussian signal 
of unit variance to the standardized forms of the arrays and and then repeating the 
statistics of Table 1. 

Prom the results of Table 2 it becomes clear that the effects of noise do not affect 
significantly the sensitivity of the SMD average in detecting the presence of the coherent 
step. 
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Table 2: Similar to Table 1 but including noise in the Monte Carlo simulations with 
signal to noise ratio - = 1.0. The full map was used in obtaining the SMD average. 



a 


Skewness 


Kurtosis 


SMD Average 


0.00 


0.01 ±0.10 


2.97 ±0.15 


0.04 ±0.31 


0.25 


0.01 ± 0.09 


2.98 ±0.13 


0.21 ± 0.34 


0.50 


0.02 ± 0.08 


2.98 ±0.14 


0.42 ± 0.39 


1.00 


0.01 ±0.10 


2.92 ±0.14 


0.97 ±0.53 



A simple way to further improve the sensitivity of the SMD statistical variable is to 
ignore a number / of boundary pixels of the SMD array, thus constructing its average using 
the Sample Mean Differences of pixels l + l, ...,n — l. From eq (24), the variance of the SMD 
for these pixels is significantly lower than the corresponding variance of the 21 pixels close 
to the boundaries. In addition, if the step is located within the central n — 21 pixels the 
SMD average may be shown to be larger than a thus further amplifying the step signature. 
For I — 150 the variance of the SMD average is reduced by about 20% (Table 3) while 
the SMD average is increased by about 20% thus allowing the detection of steps as low as 
a = 0.25 at the la level. The price to pay for this sensitivity improvement is the reduction 
of the effective pixel area where the search for steps is made. 

We have also used the SMD statistical variable for non-scale invariant power spectra 
and found that it works better for P{k) — k~'^ with < m < 1 than for m > 1. This is to 
be expected because large values of m imply larger correlations among pixels which in turn 
leads to a smaller number of effectively independent pixels and thus a larger value for the 
variance of the SMD average. 
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Table 3: Similar to Table 1 but the SMD average was obtained after ignoring 150 
pixels on each boundary of the Monte Carlo maps. The discontinuities were also excluded 
from these 300 pixels. This significantly improved the sensitivity of the SMD test. 



a 


Skewness 


Kurtosis 


SMD Average 


0.00 


0.01 ±0.10 


2.96 ±0.15 


0.01 ±0.24 


0.25 


0.01 ± 0.09 


2.95 ±0.15 


0.28 ±0.26 


0.50 


0.02 ±0.14 


2.94 ±0.18 


0.63 ±0.29 


1.00 


0.03 ±0.30 


2.78 ±0.30 


1.21 ±0.46 



It is straightforward to generalize the one dimensional Monte Carlo simulations to 
two dimensions. In that case we use the two-dimensional discrete Fourier transform as 
an approximation to an expansion to spherical harmonics. This approximation is good 
for small area maps of the celestial sphere. We used the following definition of the two 
dimensional discrete Fourier transform. 

-inn 
^ fel = l fe2=l 

refering to a n x n square lattice. In order to construct the background of scale invariant 
Gaussian fluctuations wc used g{k\^ ^2) as a complex random variable. For scale invariance, 
the amplitude of g{k\^ k-i) was obtained from a Gaussian probability distribution with 
mean and variance 

k^) = P{ki, k2) = (38) 

The corresponding phase 6^^^^^ for the (/ci, k2) mode was also determined randomly from a 
uniform probability distribution P{0ki,k2) = 2^ order to secure Gaussianity for the map 
X{i,j). To secure that the Fourier transformed map X{i,j) consists of real numbers, the 
following constraints were imposed on the spectrum g{ki, /C2) 

Img{l,l) = Img{- + 1,1)^ 
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7m^(l,| + l) = 7r„^(| + l,| + l)=0 

Tl 

g{n-i,l) = ^*(z + 2,l), ^ = o, 2 

g{hn-i) = /(l,i + 2), i^o,...,--2 
9{^:J) = /(n-i + 2,n-j + 2), i = 2,...,- + l, 

i = 2, 



The corresponding map with a superposed coherent step discontinuity was obtained 
from the standardized Gaussian map X^{i,j) as 

X'{iJ)^X^{iJ) + aj^^^^ (39) 

where 

y2-yi 



(40) 

b = yi-axi (41) 

i.e. the hne of step discontinuity j — ai + b is determined by the two random points {xi, yi) 
and {x2, ^2) of the map X{i,j). The skewness and kurtosis of the two maps are obtained in 
the usual way. For example for the standardized Gaussian map X^{i,j) we have 

1 " 

« = i^E^n^r (42) 

k = -.tx^ihoY (43) 

The SMD statistical variables is obtained by considering a set of random straight lines 
bisecting the map and for each line taking the difference of the sample means from the two 
parts of the map. For example consider a line defined by the random points {xi^yi) and 
{x2i ^2) of the map. The line equation is j = a i + 6 with a, b obtained from eqs. (40) and 
(41). The SMD obtained from this line is 

SMD = ^ - ^ (44) 
ni n2 
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where 



n n 



Si 



E E XV, j) 



(45) 



i=l j=Maxl{a i+b),l] 



n Min[(a i+h),n] 



S2 



E E x'i^,j) 



(46) 



1=1 j=i 



and Hi, 712 are the corresponding numbers of terms in the sums. For a Step+Gaussian map, 
the discontinuity the index get replaced by '. 

The average and variance of the SMD is obtained by averaging over a large number 
of random test lines (a, h) and a large number of map realizations. The results of the 
application of the three statistics (skewness, kurtosis and SMD average) on 30 x 30 scale 
invariant Gaussian maps for various values of step amplitudes a are shown in Table 4. 
Uncorrelated Gaussian noise with signal to noise ratio - = 2.0 was also included. The 
random points defining the test lines were excluded from the outermost three rows and 
columns of the maps thus reducing somewhat the variance of the SMD average. 

Table 4: A comparison of the effectiveness of the statistics considered in two 
dimensional maps. A signal to noise ratio of ^ = 2.0 was assumed in these maps. Points 
defining the line discontinuities were excluded from the three outermost rows and columns 
of the maps. 



a 


Skewness 


Kurtosis 


SMD Average 


0.00 


0.04 ±0.13 


3.00 ±0.20 


0.01 ±0.03 


0.25 


0.02 ± 0.08 


2.97 ±0.13 


0.14 ±0.09 


0.50 


0.05 ±0.14 


2.91 ± 0.24 


0.34 ±0.19 


1.00 


0.02 ±0.24 


2.95 ±0.30 


0.56 ±0.31 



The results of Table 4 are in quahtative agreement with those of Tables 1-3 and 
with the analytical results valid for the one dimensional maps. Clearly the details of the 
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one dimensional analysis are not valid in the two dimensional case and so the agreement 
can not be quantitative. The results still indicate however that the SMD statistic is 
significantly more sensitive compared to conventional statistics for the detection of coherent 
discontinuities on CMB maps. This statistic can detect coherent discontinuities with 
minimum amplitude a^m — 0.5 at the la to 2a level where a is the amplitude relative to 
the standard deviation of the underlying scale invariant Gaussian map. 



4. Conclusion-Outlook 

It is straightforward to apply the statistic analyzed in this paper to realistic data 
of ongoing experiments. Consider for example a n x n pixel sector of the COBE map 
including pixels. Let also (^)rms be the rms temperature fluctuations of the sector 
under consideration. The presence of a late long string through this sector would have 
caused a temperature step-discontinuity coherent over the whole map, with magnitude a 
given by eq. (1). If the average of the SMD over this sector is found to be very close to 
(more than la away from the SMD average value for a — 0.5), then we may conclude that 

a a 



(ST. <7^TT^<0-5 (47) 

{fjrms y-f- ) rms 

at the la confldence level, where {^)rms ^^e rms value of the purely Gaussian part of 
the fluctuations which is clearly smaller than the total {^)rms which includes the step 



discontinuity. Thus 



1 ST 

GjlVs-fsCOSe <—{ — )rms (48) 
OTT 1 



where 6 is an angle determined by the orientation of the string with respect to the observer. 
For example for {^)rms — 10~^ we obtain G/jLVs^sCOsO < 4 x 10"''' at the la level. 

Thus using the SMD statistic which is optimized to detect coherent temperature 
discontinuities on top of Gaussian temperature maps we may obtain non-trivial upper or 
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even lower bounds on the values of GjiVg^s which are highly robust and independent of 
the details of the string evolution and the resolution of the CMB maps. Application of 
this statistic on the COBE data is currently in progress (Athanasiou, Perivolaropoulos & 
Simatos 1997). 

I wish to thank G. Athanasiou, T. Tomaras and N. Simatos for interesting discussions 
and for providing helpful comments after reading the paper. This work was supported by 
the Greek General Secretariat for Research and Technology under grants UENEA 1170/95 
and 95£;A1759 and by the EEC grants CHRX - CT93 - 0340 and CHRX - CT9A - 0621. 
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Fig. 1. — A large scale coherent step- function discontinuity superposed on a one dimensional 
pixel map. 
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Fig. 2. — A standardized two dimensional pixel array of scale invariant Gaussian fluctuations. 
No step function has been superposed. 
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Gaussian+Step Scale Inv. Fluctuations. a=0.5 
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Fig. 3. — The two dimesnsional array of Figure 2 with a superposed coherent step- 
discontinuity of amphtude a — 0.5 defined by the random points (xi,yi) = (13.6,18.1) 

and (2:2,1/2) = (9.4,20.4) 



